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1. Introduction 

Simple models of intermolecular interaction have proven to be useful tools in understanding diverse 
phenomena in real fluids. Amongthese, the so-called discrete potentials that include the popular square- 
well flUlill and square-shoulder lil LtilEill potentials, or combinations of them fislEil - tiin . have received 
much attention. This is not surprising in view of their relative simplicity and versatility, which allows 
one to treat a variety of problems inclu ding , amongst others, chemical reactions liquid-liquid transi- 
tions liBjioiJiill . colloidal interactions lial , the anomalous density behaviour of water and supercooled 
fluids 134 , 35 1 and the thermodynamic and transport properties of Lennard-Jones fluids f23ll25fl . 



While the phase diagram and the thermodynamic and transport properties of the discrete-potential 
fluids have been thoroughly examined, studies of their structural properties, either theoretical or from 
simulations, are more limited. In fact, for the case where the potential is built as a combination o f sq uare 



wells and shoulders, we are only aware of the rather recent work reported in references [32, 38II39H. 

In previous papers j^lioTl, following a methodology that, although approximate, has proven success- 
ful for many other systems I41I1 . the structural properties of the square-well and the square-shoulder flu- 
ids were derived. The main aim of this paper is to use a similar methodology, referred to as the method of 
rational-function approximation (RFA), to generalize the previous results and derive the structural prop- 
erties of fluids whose molecules interact via potentials with a hard-core plus piece-wise constant sections 
of different widths and heights. 

It should be emphasized that the RFA is not a new integral equation but rather an alternative ap- 
proach to derive the structural properties of fluids in an analytic or semi-analytic way. Instead of propos- 
ing a closure relation for the Ornstein-Zernike equation, it deals directly with the radial distribution func- 
tion in Laplace space and involves some coefficients which are determined by imposing (basic) specific 
physical conditions (see references |22|,[4(l[4l|] for details). An interesting feature of the method is that, 
when applied to hard-sphere systems in odd dimensions I42H45I1 . its simplest implementation coincides 
with the exact solution of the Percus-Yevick (PY) integral equation f-iil - tiill . In the case of even dimensions, 
the mathematical problem becomes much more complicated, which is connected with the absence of a 
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simple relationship between the Laplace transform (which is a key element of the RFA approach) and the 
Fourier transform (related to the static structure factor) in even dimensions. In fact, the PY equation for 
hard-sphere fluids has no analytical solution in those dimensions. 

The paper is organized as follows. In section[2]we provide the background material required for the 
subsequent development. This is followed in section[3]by two different proposals for the RFA, both for the 
general «-step case and for the particular case of a two-step potential. Section [4] contains the results of 
our calculations for illustrative cases and their comparison with both simulation data [ 39] and with those 
that stem out of the numerical solution of the PY equation for this system. We close the paper in section[5] 
with further discussion and some concluding remarks. The consistency of the present approach with the 
results for both the square-well and square-shoulder results is proven in an appendix. 



2. Fundamental relations 



The radial distribution function g(r) of a fluid of particles interacting via a potential q>{r) is directly 
related to the probability of finding two particles separated by a distance r Isill . It can be measured from 
neutron or x-ray diffraction experiments through the static structure factor S(q), which is related to the 
Fourier transform of g{r) - 1 by 

S{q) = 1 + p f dre-^ 1 [g{r) 
G{s) - G(-s) 



1] 



\-2np 



s=iq 



where p is the number density and 



/ 



dre rs rg{r) 



(1) 



(2) 



is the Laplace transform of rg{r). The isothermal compressibility of the fluid, kj — p -1 [dp I dp] T , where 
p is the pressure and T is the temperature, is directly related to the long-wavelength limit of the structure 
factor: 

pk B TK T = S{0), (3) 

where kg is the Boltzmann constant. Thus, all the physically relevant information about the equilibrium 
state of the fluid is contained in g(r) or, equivalently, in G(s). 
Let us consider now the piece-wise constant potential 



(p{r) ■ 



£2, 



r <a, 

a < r < \\G, 
\\o < r < 



(4) 



e n , \ n -\o <r <X n a, 
0, r > A n a. 



This potential is characterized by a hard core of diameter a and n steps of (positive or negative) heights ej 
and widths (A; - A ;_i)cr (where the convention Aq = 1 is understood). Thus, A n a denotes the total range 



of q>(r). The sign of ej defines whether the y'-fh step is either a "shoulder" (ej > 0) or a "well" (ej < 0). The 
interaction potential at r = A/tr (j — 1,..., n) is repulsive if ej > ej+i and attractive if ej < ej+i (where, by 
convention, e n+ \ = 0). As usual, the density is measured by the packing fraction 77 = | pa 3 . Henceforth, 
we will take the hard-core diameter a - 1 as the length unit. 

It is convenient to define an auxiliary function F{s) directly related to the Laplace transform G{s) 
through the relation 

F{s)e~ s 



G{s) = 



l + 12r]F{s)e- 



£ (-12i7) m - 1 *[P(*)]' 



(S) 
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Laplace inversion of equation l[5) provides a useful representation of the radial distribution function: 

oo 

g{r) = r- 1 £ {-\2r 1 ) m - l f m {r-m)&{r-m), (6) 

m — l 

where f m {r) is the inverse Laplace transform of s[F(s)] m and 0(r) is Heaviside's step function. Thus, the 
knowledge of g(r) [or G[s) or S[q)] is fully equivalent to the knowledge of the auxiliary function F{s). 
The representation (6) reflects the fact that, due to the hard core at r - 1, the radial distribution function 
possesses singularities (of decreasing order) at r = 1,2,3,.... In particular, the value of g(r) at contact, 
g(l + ), is given by the asymptotic behaviour of F{s) for large 5: 



g(l + ) = /i(0) = lims z F(s). (7) 

5— 'OO 

Since g(l + ) should be finite and different from zero, we get the condition 

F{s)~s~ 2 , s^oo. (8) 
From equation (l)> it turns out that the behaviour of G{s) for small 5 determines the value of S(0): 

? 1 - S(0) 9 

G{s) = s~ z + const + — -s + 0{s z ). (9) 

2477 

Insertion of equation (|9) into the first equality of equation (5} yields the first five terms in the expansion 
of F(s) in powers of s |4oj[5si[5^1. 

It 1 ? l+2n o 2 + n A . 

F(s) = l + 5+ -s 2 + '-s 3 + -s* +&(s 5 ), (10) 

1277 { 2 12?7 24ti J 

and expresses S(0) in terms of the coefficients of s 5 and s e , namely 

24 



S(0) = — 77 3 |6 



L d 5 F(s) 




d 6 F(s) 




I 6 ds* 


s=0 


ds 6 


s=o) 



1 + 877 + 277 2 . (11) 



Equations (8) and (10) provide the behaviours for large 5 and small s, respectively, that the auxiliary 
function F{s) should necessarily satisfy. 
Let us now decompose F{s) as 

n 

F(s) = £ Rj (s)e~ W J~ 1)s , (12) 

where, as said before, Ao = 1. The aim of this decomposition is to reflect the discontinuities of g(r) at the 
points r = Aj where the potential is discontinuous. Let us denote by £;(r) the inverse Laplace transform 
of si?,-(s).Thus, 

f ? (0) = lim spRits), f'(0)= lim s[s 2 i?,(s) -^,(0)1 , j = 0,...,n. (13) 

1 5— *oo J s— *oo 

Equation (121 implies that the inverse Laplace transform of sF(s) is 

n 

;'=o 

From now on we will assume, for the sake of concreteness, that X n 2 (although the case X n > 2 can 
also be dealt with in a similar way). Therefore, equation (6) yields 

0(2-r)g(r) = r^fiir- l)8(r- 1)0(2- r) 



r ~ l L " A ;) ( r " A i)©( 2 " r )- (IS) 
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As a consequence, 



gap - gwj) = Ujm, g'ap-g'ap = ± km - Ujm 



j = 0,...,n. 



(16) 



On the other hand, the cavity function y(r) = g(r)e' p ^ lkBT and its first derivative y'(r) are continuous at 

7 = 1,..., ra, 



r = Xj 15711 . This means that 



g(AT)e^ = g(At)e^«, g'(AJ)e^ = g'U})e^ 
or, according to equation (15), 

E 

;"=0 



^■(0) = \efcr<J+r) - 1] £ ^(A;-A,), ; = 1, 



f;.(0) = [e^- e ^-l]2:^(A;-A i ), ; = 1 

i=0 



(17) 



(18) 



(19) 



In the low-density limit, one has g(r) -» e ^ (r) , so that 



limF(s) = s'Vlim G(s) 



r/— 



77— 



s -3 £ e - (A J- 1)s (1 + Ajs) (e-*/ +1 - e"^) , 



(20) 



7=o 



where, apart from e„+i = and Ao = h it is understood that eo = 00. Comparison with equation (12) yields 
limJ?7(s) = s _3 (l-i-A7s)(e^ +1 -e^j, j = 0,...,n. (21) 

3. Construction of the rational-function approximation 

3.1. The general n-step case 

According to the scheme presented in the preceding section, the full knowledge of the radial distribu- 
tion function g(r) for the rc-step potential (4) amounts to a prescription for the functions Rj(s), such that 
the function F{s) obtained through equation (12) fulfills the conditions (8) and (To) . 

Now we assume the following rational-function approximation (RFA) for Rj (s): 



Rj(s). 



1 



Aj + Bjs 



1277 1 + S 1 s + S 2 s z + S 3 s : 



3 ' 



j = 0,...,n. 



(22) 



Since the degree difference between the numerator and denominator of Rq{s) is equal to 2, the form (22) 
ensures the consistency with equation (8). In fact, according to equation (13) , 



1 Bi , l 
12?7 S3 1 \2r\ 



'Aj_ BjS 2 \ 

I S 3 



s 2 



7 = 0,..., n. 



Insertion into equation (16) yields 



7=0,..., n, 
(4/ *7«2 B7 \ 



S3 S3 S3 Ay ; 



7 = 0,... ,n. 



(23) 



(24) 



(25) 
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For r > 0, application of the residue theorem provides the inverse Laplace transform of sRj (s) as 

f (r) = -— £ , s a e'° r , (26) 

S; 1217^1 Si +2S 2 s«+3S 3 4 

where s a (a = 1,2,3) are the three roots of the cubic equation 1 + Sis a + S2S 2 a + S3S 3 = 0. 

The approximation (22) contains 2{n + 1) + 3 = 2rc + 5 parameters to be determined. The exact expan- 
sion (10) imposes five constraints among them: 

n 

A = l-Y J A } , (27) 

i=i 

Si=B -(l + C (1) ), (28) 
S 2 = ^(1+2C (1) + C (2) )-B , (29) 

S 3 = ifi - — - - (1 + 3C a) + 3C (2) + C (3) ) , (30) 
2 12tj 6 l ; 

(1 + 2r))B = 1 + C (1) + I (l + 4C (1) + 6C (2) + 4C (3) + C (4) ) , (31) 



where 

C m = £ UyU; - l) fc - JfcflyCA,- - l) fc_1 | . (32) 

Thus, the five coefficients ^4o, Bo> Si, S2 and S3 are given by equations (27)-(32) as linear combinations 
of Ai and Bj (j = l,...,n). In the case of a hard-sphere system (formally, n — 0) the problem is closed 
and one recovers ^4ll, [55L [HSl the exact solution of the PY integral equation I46T - I481] . On the other hand, 2n 
parameters remain to be determined ifn^X. 

Before dealing with this problem, let us see that the proposal (22) is consistent with the exact form 1 121) 
in the zero-density limit. This implies that one should have 

limSi=0, lim77S 2 = 0, lim77S 3 = — -, (33) 
limA = e"^ ei , (34) 

r/— 

limflo = e _/3<ri , (35) 
7;— 

limAf = e~^ +1 -e~^, i = l,...,n, (36) 
;j-o ' 



limBj = Xj (e"^ +1 - e"^j , j = l,...,n. 



(37) 



First, note that, if equation (36) is satisfied, then equation (27) reduces to equation (34) . Next, equa- 
tions (36) and (37) imply that lim^o C' 1 ' = Aq - 1. Thus, equation (31) reduces to equation (35) . Finally, 
taking the limit 77 — ► in equations (28)-(30) one obtains equation (33) . This proves that the approximate 
form (22) becomes exact in the limit 77 — ► 0, provided that the 2n coefficients Aj and Bj (j = 1, . . . , n) are 
consistent with equations (36) and (37) . 
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3.1.1. RFA1 

As a first proposal to determine the 2n coefficients Aj and Bj with j = l,...,n, we may discard the 
density dependence of Aj (j = 1, . . . , n). As a consequence, those coefficients are fixed at their zero-density 
values, namely 

Aj=e~ Pe } +l -e~ Pe 1, j = l,...,n. (38) 

Insertion of equation (38) into equation (27) implies that Aq is also given by its zero-density expression, 
equation (34). The coefficients Bj (j - \,...,n) are determined from the n constraints (18) stemming from 
the continuity of y(r). Making use of equations (23) and (26) yields 



?L = Ubi-'j^ - 1] £ n S f X}S \ c , Z {A t + B iSa )e-^ s ", j = l,...,n. (39) 



These conditions should be enforced numerically. We will refer to this variant as the RFA1 method. In 
the one-step case (n = 1), it coincides with the one previously proposed for square-well fllL I40I [4ll [s^jsoll 
and square-shoulder 12211 fluids. The RFA1 reduces to the PY solution for hard spheres and sticky hard 
spheres in the appropriate limits fzzL [4oh . 



3.1.2. RFA2 

A more sophisticated version (here referred to as RFA2) consists in replacing the simple prescrip- 
tion (38) by the enforcement of the continuity of y'(r) via equation (19) . Again, taking into account equa- 
tions (23) and (26) , one gets 

^_^ = [^^i)_i]f f l]S " - ZiA i + B t s a )e-^, j = l,...,n. (40) 

Now the problem requires the numerical solution of the 2n coupled transcendental equations (39) and 
(40) . instead of the n equations (39) required in the RFA1. 

While the RFA2 is internally more consistent than the RFA1, it is not necessarily more accurate. As we 
will illustrate in section[4| the requirement of a rather subtle continuity condition of y'ir) at r = Aj may 
force a radial distribution function with features less realistic than those found. 

In either case, once the Laplace transform G[s) is fully determined, its numerical Laplace inver- 
sion (ill] yields g(r). 



3.2. Particularization to two-step potentials 

Let us now adapt the previous scheme to the case n — 2. Then, 

F[s) ^ R {s) + R 1 (s)e~ Ul ~ 1)s + R 2 (s)e~ U2 ~ 1)s , (41) 



where the three functions Ro{s), R\{s) and Rzis) have the rational-function form (22) . There are 9 coef- 
ficients to be determined. In the simpler RFA1 approach, Aj are fixed by their zero-density values, i.e., 
according to equation (38) . 



A = e-P £l , A 1 = e~P e2 -e-f kl , A 2 = \-e-P e2 . (42) 

Next, Si, S2 and S3 are given as linear combinations of By by equations (28)-(30), where 

C (fc) = A x CAi - 1) k + A 2 (A 2 - 1) k - kBi (Ai - 1) k ~ l - kB 2 (A 2 - 1) fc_1 . (43) 

Moreover, equation (31) expresses Bo in terms B\ and B 2 . Thus, only two equations are needed. They are 
given by equation (39) , i.e., 

J. = L^-^-l Z " 2 (A + B s a ), (44) 

S3 L J a= iSi+2S 2 s a + 3S 3 5^ 



23602-6 



Piece-wise constant potentials 



|i = (e^ 2 -l) £ S ^ll^——\ {Ao+BoSa)e -^ + {A +B s )e -hs a } 

These two transcendental equations close the problem. 

In a more sophisticated RFA2 version of the approximation, equation {42) is replaced by 



S3 Sj 1 £? 1 Si+2S 2 S a + 3S 3 si 



A 2 B 2 S2^ Uez _ A ^ s 2 a e^ 

;2 



(45) 



Ao = l-Ai-A 2 , (46) 

and equation {40), i.e., 



1 E 7 — ^-Uo + Bo%), (47) 



f 6 ^ 2-1 ) E g ^/" 6 2 [(^ + Jb^)e-*» + (Ai+Bi5 8 )e- A ^]. (48) 

S3 S^ v a— 1 Si + 2S2 Sa + 3S3 J 

These two equations should be solved in conjunction with equations {44) and {45) . 

It is shown in the appendix that the RFA for the two-step potential reduces to the one for the one-step 
potential in the limit e\ = e 2 as well as in the limit ei — > 00. 



4. Illustration of the approximation for two-step potentials 

In this section we illustrate our proposals RFA1 and RFA2 in two representative examples of two-step 
potentials. First, we consider a square-shoulder plus square-well (SS+SW) potential (i.e., e\ > and e 2 < 0) 
withei = \e 2 \, Ai = 1.25 andA 2 = 1.5. The second case corresponds to a shifted square-well (sSW) potential 
(i.e., ei = and e 2 < 0) with Ai = 1.25 and X 2 - 1.5. Both potentials have recently been analyzed by 
Barcenas et al. fiiil by means of exchange replica Monte Carlo (MC) simulations, with special emphasis on 
the coexistence curves between the vapour and the condensed phase (liquid or solid) and the structural 
properties of the condensed phase at coexistence. 

The chosen state points are per 3 = 0.421, fc B T7|e 2 | = 0.64 for the SS+SW system and per 3 = 0.427, 
k^TI\e 2 \ = 0.76 for the sSW system, both lying on the condensed branch of the respective coexistence 
curve 13911 . The radial distribution function g(r) at a subcritical temperature and at the density of the 
coexisting liquid represents a rather stringent test of the theoretical approach developed in this paper. 
To put the results in a proper context, we have also numerically solved the PY integral equation by an 
iterative method 12211 . The comparison between RFA and PY is especially relevant since both theories 
coincide in the hard-sphere and sticky-hard-sphere limits I22II40I1 . 

Figures[T]and[2]show the radial distribution function and the cavity function for the SS+SW and sSW 
systems, respectively, as obtained from MC simulations liill . the PY numerical solution and our RFA1 and 
RFA2 approaches. We observe from the figures that the RFA1 and RFA2 curves are practically indistin- 
guishable in the region r > A 2 , where they describe the main oscillating trends of the true distribution, at 
least in a semi-quantitative way. Inside the well (Ai «r« A 2 ), however, both versions of the RFA differ. In 
the SS+SW case (figure [1), the RFA2 prediction is excellent, visibly differing from the MC data only near 
the inner radius of the well, r -X\, while the RFA1 exhibits an accused convex shape; as a consequence, 
in the shoulder region (1 =S r =£ Ai) the RFA1 predicts a g{r) with positive slope. In the sSW case (figure[2), 
on the other hand, it is the RFA2 that presents a convex shape inside the well (Ai =£ r =£ A 2 ) and hence 
exhibits a positive slope in the shift region (1 s£ r < By contrast, the RFA1, which (as clearly seen in 
the bottom panel of figure[2) has a discontinuous first derivative at r = A 2 , presents an almost linear form 
in the regions 1 ^ r ^ X\ and X\ ^ r ^ A 2 with a qualitatively correct behaviour. Note, however, that the 
RFA2 is very accurate in the region inside the well and near the outer radius (1.4 «r« A 2 ). 

As for the PY theory, it is obvious from figures [1] and [2] that it performs worse than any of the RFA 
results for r^X\ and than either RFA2 (SS+SW system) or RFA1 (sSW system) for 1 ^ r =£ X\. This remark- 
able result (given that the PY solution requires full numerical work, in contrast to the semi-analytical 
character of the RFA) was already observed for square-shoulder fluids I22I1 . 
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Figure 1. (Color online) Radial distribution func- 
tion g(r) (top panel) and cavity function y(r) (bot- 
tom panel) for a square-shoulder plus square-well 
potential (e\ > and e 2 < 0) with ei = fal, Ai = 
1.25 and A 2 = 1.5 at per 3 = 0.421 and knTI\e 2 \ = 
0.64. The thick solid lines are MC results (Hi, while 
the dashed-dotted, dashed and thin solid lines cor- 
respond to the PY, RFA1 and RFA2 theories, respec- 
tively. 



Figure 2. (Color online) Radial distribution func- 
tion g(r) (top panel) and cavity function y(r) (bot- 
tom panel) for a shifted square-well potential (ei = 
and £2 < 0) with Ai = 1.25 and A2 = 1.5 at per 3 = 
0.427 and knTI |e 2 1 = 0.76. The thick solid lines are 
MC results 13911 . while the dashed-dotted, dashed 
and thin solid lines correspond to the PY, RFA1 and 
RFA2 theories, respectively. 



5. Conclusions 

In summary, in this paper we have proposed a method of deriving the structural properties of a 
particular kind of discrete-potential fluids, namely the ones in which molecules interact via a potential 
consisting of a hard-core plus rc-step constant sections of different heights and widths. The method is 
based on assuming rational-function forms for n + 1 functions Rj [s) defined in Laplace space, the coeffi- 
cients being constrained by physical consistency conditions. Two approximations were considered: one 
in which some unknown constants are fixed at their zero density value (RFA1) and one which enforces the 
continuity of the derivative of the cavity function (RFA2). In the former case, one has in the end to solve 
(numerically) n coupled transcendental equations while in the latter, which is internally more consistent, 
the set to be solved is made of 2n coupled transcendental equations. 

The illustrative cases of the SS+SW and sSW potentials that we examined indicate that our approach, 
being semi-analytical in nature, is a reasonable compromise between simplicity and accuracy. For the 
SS+SW potential, the RFA2 approximation gives the best results, whereas the RFA1 is superior in the 
sSW potential. On the other hand, in these illustrative cases the present development also represents a 
clear improvement over the solutions of the corresponding PY integral equations, which require harder 
numerical labor. 

These further examples of the adequacy of the RFA method as compared with the results of the PY 
equation reinforce the notion that such methodology is a valuable alternative to the integral equation ap- 
proach for the derivation of the structural properties of fluid systems. This conclusion is of course based 
on a limited analysis. One could easily argue that the comparison with the results stemming out of the 
Ornstein-Zernike equation with either the Rogers-Young or the hypernetted chain closures would be a 
more solid test ground. However, these approximations require much harder numerical work, while our 
approach yields the explicit s-dependence of the Laplace transform G(s). In any case, we find it interesting 
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to compare the performance of the present development with the results of some other simple approxi- 
mations such as the hybrid mean spherical approximation considered in reference 1 3 211 or the first-order 
mean spherical approximation used in reference I38I1 . We plan to undertake such comparisons in our 
future work. 
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A. Consistency with the one-step (square-well or square-shoulder) case 

Let us denote by g(r;ei/f: 2 , Ai, A 2 , T*,rf) the radial distribution function in the case of a two-step po- 
tential with reduced temperature T* = /cb 77e 2 and packing fraction 77 = | p, where the hard-core distance 
is a = 1. Analogously, we denote by g(r; A, T*,T]) the radial distribution function of a one-step (square-well 
or square-shoulder) potential of unit hard-core distance and relative width A with reduced temperature 
T* (in units of the step height) and packing fraction 77. 

The two-step problem should reduce to the one-step problem in certain limits. The aim of this ap- 
pendix is to prove that the RFA2 method fulfills these consistency conditions. A similar and simpler proof 
can be worked out for the RFA1 method. 



A.1. Case ei = e 2 

An obvious physical condition is that the two-step system becomes equivalent to a one-step system in 
the limit e\ — ► e 2 , i.e., 

g(r; LAi.Az, T*,T]) = g(r; A 2 , T*,1]). (49) 

In Laplace space, 

G(s;l,Ai,A 2 ,r*,ij) = G{s;A 2 ,T*,T]), (50) 

F(s;l,Ai,A2,r*,r?) =F(s;A 2 ,r*,7?). (51) 
Taking into account equations (12) and (22) . one should have 

A (l,Ai,A 2 , T* ,77) = ,4o(A 2 , T*,T]), (52) 

Ai(l,Ai,A 2 ,T*,?7) =0, (53) 

A 2 a, Ai, A 2 , T* ,77) = Ai (A 2 , T* , T]), (54) 

£ (l,Ai,A 2 ,r,77) = .Bo(A 2 ,r\77), (55) 

£i(l,Ai,A 2 ,r*,77) = 0, (56) 

B 2 (l,Ai,A 2 ,r*,7 7 )=B 1 (A 2 ,r*,T7), (57) 

S 1 (l,Ai,A 2 , r*,??) = Si(A 2 , r* ,77), (58) 

S 2 (l,Ai,A 2 , T*,r]) = S2U2, T*,T]), (59) 

S 3 (oo,Ai,A 2 , T*,r]) = S 3 (A 2 , r*,77). (60) 

Equations f52l-j54l are in full agreement with equation (27) . It remains to prove that the above equations 
are consistent with equations (28VJ3TV d44> , (45) . (47) and (48). 

For simplicity, henceforth we will omit the arguments of the barred and unbarred quantities. From 
equation (43) , 

c (fc) _ £{k) _ ^ i(A2 _ 1} fc _ k B^X 2 - l) k ~\ (61) 
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where use has been made of equations (53")-(57). This ensures that equations d28t— d3Tt are satisfied by 
the one-step system, provided they are satisfied by the two-step system in the limit ei — » ez. Next, the 
transcendental equations (44) and (47) are trivially satisfied by equations (53) and (56). Finally, the tran- 
scendental equations (45) and (48) become the same for the two- and one-step cases. 



A.2. Case ei — oo 

In the limit e\ — <■ oo the two-step potential becomes equivalent to a one-step potential of hard-core 
distance X\ and relative width A = A2/A1 with a packing fraction fj = ^A 3 , i.e., 

g(r;oo,Ai,A 2 .r*,77) = g{rlX\; A, T*,fj). 

In Laplace space, equation (62) translates into 

G{s; 00, Ai , A 2 , T* , 77) = X\ G(Ai s;X, T*,fj). 

From equation (5), 



F(5;oo, Ai, A2, r*,?7) = A 3 e 
According to equations (12) and (22), one should have 



3 -(Ai-l)s 



F(Ais;A, T*,fj). 



Aq (00, A 
A\ (00, A 
^2 (00, A 
_Bo (00, A 
B\ (00, A 
B2 (00, A 

51 (00, A 

52 (00, A 

53 (00, A 



,A 2> r,77) = u, 
,A 2 ,r*,?7) = ^o(A ( T*,fi), 
,Az,T*,Tj)=A 1 tf., r,jj), 
,A 2 ,r*,?7) = 0, 
,A 2 ,T*,j 7 )=AiBo(A ) r*,j7), 

,X z ,T*,rj)=\ l B 1 {X,T*,fi), 
,A 2 ,r*,7 7 )=AiSi(A, r,7)), 

,A 2 ,r*,7 7 )=A?S 2 (A, r,?7), 

,A 2 ,r,7 7 ) = A?s 3 (A, r,77). 



(62) 

(63) 
(64) 

(65) 
(66) 
(67) 
(68) 
(69) 
(70) 
(71) 
(72) 
(73) 



Again, equations (65)-(67) are consistent with equation (27) . Let us now check that equations (65)-(73) 
are consistent with equations (28)-(3l), (44) . (45) , (47) and (48) . 

As before, for simplicity we omit the arguments of the barred and unbarred quantities. Taking into 
account the definition (43) , one has 



C [k] = Ai(A-l) k - kBi[X-l) k ~ 1 . 
According to equations (65)- (70). we get that, in the limit ei -» 00, 

C m = (Ai - l) fc + Aj [(A 2 - l) k - (Ai - l) fc ] - kB Xi{Xi - l) fc_1 - fcBiAi(A 2 - l) fc "\ 
where we have made use of the property Aq + A\ — 1. From equations (74) and (75) one has 

l + C m =Xi{l + C m -B ), 

1 + 2C (1) + C (2) = X\ (l + 2C (1) + C (2) - 2B ) . 

1 + 3C™ + 3C (2) + C (3) = Aj (l + 3C (1) + 3C (2) + C (3) - 3B ) , 

1 + 4C (1) + 6C t2) + 4C (3) + C (4) = At (1 + 4C (1) + 6C (2) + 4C (3) + C (4) - 4B ) . 

Equations (76)-(78), together with equations (28)-(30). prove the consistency of equations (7T)-(73). _ 
spectively. Next, equations (76) and (79) prove that, if (31) is satisfied by the unbarred parameters, so it is 
by the barred parameters. 



(74) 
(75) 

(76) 

(77) 
(78) 
(79) 
re- 
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Now we consider the transcendental equations that need to be additionally solved. In the case of the 
one-step potential, equations (39) and (40) reduce to 



B\ I tit' \ ^ s a e {X - 1)s " 
S 3 v Si+2S 2 s a +3S 3 s^ 

Ai BiS 2 i jit* £e a ~ l)s « 

03 o 3 v ' a=l Si +£02S a +i$-iS a 

respectively, where s a (a = 1,2,3) are the three roots of the cubic equation 1 + S\s a + S 2 Sa + S3S„ - 0. In 
the case of the two-step potential, multiplication of both sides of equations (44) and (47) by e~P Cl shows 
that they are identically satisfied in the limit ei — 00 if Aq = Bo = [cf. equations (65) and (68)]. As for 
equations (45) and (48), they become 

T =e UTT^ -rz-zlAi+BiSa), (82) 

S3 V ^ 1 Si+2S 2 5 a + 3S 3 4 



A 2 B 2 S 2 ( ut* 1^ s 2 (A 2 -Ai)s„ 

— --^r L = e 1/r -1 }Y 5-(^i + BisJ. (83) 

S3 Sj y )^ l S 1 +2S 2 s a + 3S 3 s z a K 

Taking into account the equations (66) . (67) and (69")-(73). as well as the fact that s a - X\ s a , it is straight- 
forward to check that equations (80) and (81) are equivalent to equations (82) and (83) . respectively. This 
closes the proof of equation (62) . 
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Ha6/ii/i>i<eHHfl pai^ioHa/ibHMMi/i tyyHKU,\nMV\ pjin ruii/miB, mo 
B3aeM0fliK)Tb Hepe3 BiflTHHKOBO-CTa/ii noTemjia/ii/i 

A. Camo^, C.B. rOcnP, M. /lonec f\e Tapc@ 

1 ct>i3H4Hnw c|>aKy.fibTeT, YHiBepcuTeT EKCTpeMaflypn, E-06071 BaflaRoc, IcnaHin 

2 Llemp eHepreTH4Hwx ,qoc/ii,q>KeHb, HauiOHa/ibHui/i aBTOHOMHuti yHiBepcmeT MeKcnKW (U.N.A.M.), 
Mopeyiboc 62580, MeKCHKa 

BHK0pHCT0ByKJ4H (HaniBaHa^iTH4HWii) MeTOfl Ha6/lW>KeHHfl pauiOHa/lbHUMW (fiyHKL^mMW, OTpniviaHO CTpyKTypHi 

B/iacTHBOCTi n/inHiB, MO/ieKy/in akhx B3aeMOfliK)Tb 4epe3 noTeHL(ia.nn i3 TBepflOio cepueBUHOK) Ta n cxoflUHKa- 
mm pi3HOi LunpwHW Ta bucotw. Pe3y/ibTa™ i/iiocTpyioTbCfl Ha npuK/iaflax noTei-mia/iy npaMOKyTHnw ycTyn n/iroc 
npHMOKyTHa HMa i noTeHL(ia/iy 3MimeHa npflMOKyTHa flMa, i nopiBHKJKnbca 3 pe3y/ibTaTaMH cv\wym\\\\A, a Tai<o>K 
3 pe3y/ibTaTaiviw, mo OTpnMyiOTbCfl 3 (4hoiobwx) po3B'fl3KiB iHTerpa^bHoro piBHAHHfl nepKyca-CBiKa. 

K/iK>MOBi oiOBa: fli/iCKpeTHi noTeHU,iajw, MOflejib npuMOKyTHoro ycryny, MO/s,enb npnMOKyrHo'i rmm, paflia/ibHa 
<pyHKU,iR po3nofliny, (pyHKLfin nopo)KmiHi/i, Ha6mvKeHHR pau,iOHajibHi/iMti (pyHKLfinMH 
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